k = 8;
z = [-1 + 1j, -1 - 1j];
p = [0, 0, -5, -6, 1j, -1j];
Gs = zpk(z, p, k)